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ABSTRACT 


Two  u«  methods  were  developed  for  predicting  projeceilo  shape 
which  yield  minimum  total  drag  at  supersonic  speeds.  The  first  tech** 
nlque  is  an  Eulerian  scheme  that  uses  modified  Newtonian  theory  end 
Prandtl-Meyer  expansion  for  pressure  drag  with  Van  Driest  skin  friction 
and  semi-empirical  base  drag  prediction.  The  second  scheme  Iterates 
body  coordinates  with  the  second-order  shock-expansion  theory  and  the 
same  skin  friction  and  base  drag  methods  to  minimize  the  total  drag. 

A  different  shape  is  determined  for  each  length-to-diameter  ratio  and 
Mach  number. 

The  first  technique  was  found  to  calculate  a  reasonably  accurate 
optimal  shape,  but  did  not  predict  accurate  drag  coefficients.  It  was 
found  that  the  modified  Newtonian  theory  plus  Prandtl-Meyer  expansion 
predicted  pressure  drag  coefficients  much  too  low  whereas  the  second- 
order  shock-expansion  method  gave  good  results.  The  second  technique 
predicted  both  accurate  optimal  shapes  and  drag  coefficients.  Optimal 
shapes  were  predicted  using  the  second  technique  for  Mach  numbers  2-5 
and  length-to-diameter  ratios  of  4,  5,  and  6.  They  were  found  to  com¬ 
pare  well  with  experimental  data. 


ACKNOWLEDGEMENT 


The  author  would  ilk*  to  axpreaa  hie  appreciation  for  tha  guidance 
and  help  given  by  Dr.  Fred  R.  DeJarnatta  of  North  Carolina  State  Unlver- 
aity.  Hia  help  in  developing  tha  Eulerian  Optimization  waa  indeapenaible 
and  hia  encouragement  waa  continual.  Alao  thanka  are  given  to  Dr.  Frank 
G.  Moore  and  Dr.  Charlea  M.  Blackmon  of  the  Naval  Surface  Weapona  Center 
for  their  support  during  this  effort. 


v 


CONTENTS 


IMS. 

INTRODUCTION . 1 

FIRST  PREDICTION  METHOD . .  .  .  3 

EULER IAN  OPTIMIZATION  .  7 

I 

Forebody . 7 

Afterbody . II 

Eulerlan  Results . 17 

OPTIMIZATION  USING  SECOND-ORDER  SHOCK- EXPANSION  PREDICTION.  ...  19 

Eulerlan  Technique . 19 

Second-Order  Shock-Expansion  Prediction  Method  .......  19 

New  Optimization  Method . 22 

Results . 26 

CONCLUSIONS  AND  RECOMMENDATIONS  .  28 

REFERENCES . 29 

GLOSSARY  OF  TERMS . 40 

DISTRIBUTION . 42 


vi 


FIGURES 


Number 

1. 

2. 

3. 

4. 

5. 

6. 

7. 

8. 

9. 

10. 


Title  Page 

Mean  Bass  Pressure  Coefficient  30 

Projectile  Coordinates 

Comparison  of  Drag  Coefficient  vs  L/D  for 

Euler ian  Optimisation  and  Reference  5  32 


Optimum  Shape  Comparison,  Eulerian  Optimization 
and  Reference  5,  Mach  No. *3,  L/D«5 
Drag  Coefficient  V3  x^x^  fcr  h/P*4, 

New  Optimization 

Drag  Coefficient  vs  x^x^  for  L/D-5, 

New  Optimization 

Drag  Coefficient  vs  x£/xf  for  L/D-d, 

New  Optimization 

Optimal  Body  Shape  for  Mach  No, *3,  x^/x^-0.7, 
L/D»5,  New  Optimization  and  Reference  5 
Comparison  of  Drag  Coefficient  vs  x^x^  f°r  Mach 
No. »3,  L/D»5,  New  Optimization  and  Reference  5 


33 

34 

35 

36 

37 

38 


Drag  Coefficient  vs  s  /x.  for  Mach  No. -3, 

c  t 

New  Optimisation 


vit 


INTRODUCTION 


The  need  for  optimal  body  design  is  minimizing  total  drag  has 
bean  generated  by  the  Navy's  requirement  for  prcjectilea  to  have  longer 
range,  shorter  flight  times,  and  higher  terminal  velocity.  The  work  of 
Moore1  in  1959  analyzed  optimal  projectile  shape  using  mostly  empirical 
techniques.  Results  of  this  study  Indicated  that  range  of  current  projec¬ 
tiles  could  be  increased  by  more  than  fifty  percent  using  aerodynamic 
design  considerations.  Because  of  high  experimental  costs,  it  is  desired 
that  optimal  body  shapes  be  generated  by  cheaper  analytic  means.  T  ree 
main  contributes  of  drag  must  be  predicted  in  order  to  evaluate  total 
body  drag.  Th«y  are  pressure  or  wave  drag,  skin  friction, and  base 
drag.  The  major  portion  of  analytical  *-ork  has  been  in  the  prediction  of 

optimal  forebody  shapes  by  minimizing  pressure  drag.  Minimum  wave  drag 

2  3 

shapes  were  found  by  von  Karmen  ,  using  slender  body  theory.  Cole  using 

4 

Newtonian  theory,  and  Miele  ,  who  included  skin  friction  drag  with  pres¬ 
sure  drag  calculations.  The  above  optimization  studies  have  led  to 
configurations  which  adequately  predicted  optimum  supersonic  nose  shapes 
but  neglected  base  drag  contributions.  The  work  of  Hager,  et.  al.^ 
attempted  to  define  optimal  projectile  shape  including  total  drag 
analytically.  However,  lAien  compared  to  experiment,  the  drag  predicted 
was  found  to  be  low. 

The  ebjeo:  of  this  effort  waa  to  create  e  more  accurate  technique  of 
analytically  predicting  minimized  total  drag  body  shapes.  The  super¬ 
sonic  regime  (Mach  numbers  2-6)  was  chosen  since  projectiles  were  the 
bodies  of  interest.  The  approach  was  to  first  try  a  different  optimize  x 


schaaa  Chao  that  of- Kafaranca  5  vhila  a  till  using  tha  um  drag  pradiction 
nathodology-  K  aora  accurata  praaaura  drag  pradiction  tachnlqua  vaa  than 
triad  to  furthar  improva  optimization.  K  third  optimization  achama  vaa 
finally  found  that  java  aora  accurata  raaulta  although  it  vaa  conpu- 
tationally  aora  tlaa  consuming. 


FIRST  PREDICTION  METHOD 


An  Kulerian  optimisation  schema  vu  first  triad  on  ths  drag 
predictive  tschnlquss  of  Rafernce  5.  Ths  Eulerlan  optimization  is 
similar  to  that  of  Mlala  in  Reference  4.  Ths  drag  prediction  mathodologj 
of  Rafernce  S  usas  aodlflad  Newtonian  prassura  distribution  plus 
Prand tl -Mayor  expansion,  Van  Drlast  turbulant  skin-friction  analysis, 
and  a  seml-eapirlcal  basa  drag  prediction.  Tha  total  drag  coefficient 
Is  defined  by 


s  -  *  t  vx>  ctr' <x»  “ + cf.  sf  -  cx,  ( i  f  (i> 


where 


Cp(x)  is  the  pressure  coefficient  predicted  by  isodifled  Newtonian  theory 


plus  Prandtl-Meyer  expansion,  C_  is  tha  flat-plata  turbulant  skin- 

friction  coefficient,  and  C  is  tha  basa  pressure  coefficient.  Tha 

PB 

forebody  pressure  coefficient  is  found  using  modified  Newtonian  theory 


C  ■  C  sin*  8  (2) 

p  po 


where 


C  is  the  stagnation  pressure  coefficient  behind  a  normal  shock 
po 

defined  by 


<Y  +  1)  *& 
2 


T»1 

2yM»  -  (r-D 


1 

Y-l 

-  1  (3) 


and  9  is  the  body  slope  with 


3 


T*  (*)  -  UB  # 


(4) 


The  stagnation  p  rata  but*  calculation  1  limited  Iron  a  bluntad 

ooM  co  tha  point  of  maximal  thickness.  At  tha  point  of  maximum 

thlcknaaa  Equation  (2)  glvaa  C  -  0  loading  to  small  lnac- 

P 

curias. 

tha  afterbody  prssaura  calculation  Is  calculated  from  tha 
Fraud tl- Meyer  expansion 


dp  ypM2 
d9  m2-i 


(5) 


This  expression  Is  limited  to  negative  slopes  of  less  than  8°  on  the 
afterbody. 

The  skin-friction  prediction  is  from  Reference  6  which  assumes  a 
fully  turbulent  boundary  layer.  The  mean  skin  friction  coefficient 
for  a  flat  plate,  Cf  Is  found  through  iteration  of  the  following 
equation 


4 


Tha  variable  a  in  Equation  (6)  is  ths  povar  in  the  power  viscosity 


lav 


M 


a 


(7) 


and  for  air  a  Is  0.76.  A  Prsndtl  number  of  unity  and  a  zero  pressure 
gradient  in  tbs  fully  developed  turbulent  boundary  layer  are  assumed  In 
the  relations  above.  The  freestreaa  Reynolds  nunher  is 


Re 


*1 


y„ 


(8) 


where 

£  is  the  total  configuration  length.  The  temperature  ratio  TV/T^ 

Is  found  assuming  an  adiabatic  wall  and  by  Introducing  the  turbulent 
recovery  factor 

Tw  , 

T"1  +  RT2«'  (' 

OK 


The  turbulent  recovery  factor  varies  approximately  as  the  cube  root 
of  the  Prandtl  maker  i.  that 


*T  -  3iJPt  .  (10) 

To  compensate  for  the  assumption  of  Pr«i  In  the  Vac  Driest  method, 
the  actual  Prandtl  number  of  air,  Pr*»0.73  was  used  In  Equation  (10). 
Thus  Equation  (9)  becomes 

Tv  T  1  2 

Y  -  1  +  0.9  -j-  tf  . 
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(11) 


One  can  now  combine  Equation*  (8)  and  (11)  with  (6)  to  solve  for  C„  . 

The  Newton-Raphson  method  Is  used  to  calculate  iq  Equation  (6). 

The  base  drag  is  calculated  using  a  semi-eaperical  technique 
developed  by  Moore^.  A  mean  curve  of  ^rperimental  base  pressure  data  is 
given  In  Figure  1.  This  data  assumes  a  fully  developed  turbulent 
boundary  layer  ahead  of  a  long  cylindrical  afterbody.  The  effect  of  a 
boa tail  significantly  alters  base  pressure  and  oust  be  accounted  for. 

The  empirical  equation  used  is 


fir.-c  fi*-]3  . 

I  dr  J  PBA  l  dr 


c„  -  -c 
db  PI 


Equation  (12)  can  bo  used  for  all  Mach  numbers  where  C  Is  the  base 

PBA 

pressure  given  in  Figure  1. 


6 
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The  type  of  body  being  optimized  Is  of  the  general  configuration  shown 
in  figure  2.  It  consists  of  a  b inured  nose,  and  boatail  with  a  dis¬ 
continuity  9-t  the  corner.  A  mtrlmuw  body  radius  and  length  to  diameter 
ratio  (L/D)  are  constraints.  Using  (14) f  (13),  and  (16)  consider  the 
first  variation  of  equation  (13)  as 

x. 


C_r2 
D  max 


\  , 
/ 


'  *i 

5  rdx  + 


Fr  -  ~  Fr’ 
dx 


6  rdx  + 


c+ 


[£F  -  r?5r,)  Sx  +  Frl6^j^  -  pF-r‘Fr')6x  +  AF/Sr] 


CF°  ?1  5ri  “  =  C?AB  r!  irf 


(17) 


where 

the  subscripts  C-  and  C+  denote  conditions  immediately  before  and  after 
the  corner  at  X£  and  5r  represents  the  distance  between  the  external 
and  comparison  arc  as  defined  in  Reference  4.  For  minimum  drag,  the 
variation  of  drag  in  equation  (17)  equals  zero  so  now  consider  the 
right  side  of  the  equation  as  distinct  parts. 

For  the  integrand  of  the  integrals  in  Equation  (17)  to  be  zero, 


F 

r 


_d_ 

dx 


(Fr’)  -  0 


but  since  F  is  a  function  of  r  and  r' , 


then  from  Miele 


4 


(18) 


F  -  r’Fr*  •  constant  ■  -C  . 

u 

Recalling  chat 


F  - 


+  rfyjr'2  +  l 


f 


(19) 
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C  -  C  (r’> 
P  P 


and 


V 


d(C  r’>  C,r’ 

*  — tBi —  +  -  f 


dr' 


Vr* 2  + 1 


Then, 


F  -  r-yr.  -  j^.  +  Cf  ^-Tt]  r  -  r 


dCP 

CP  +  *'  d^ 


,  c, 

rr  f 
/r '  z  +  1 


SP^TT 


,  dCP 
-  r»2  _ L 

dr* 


and  further 


(20) 


dC_ 

■  •2  _J? 

dr’ 


/r'*  +  1 


constant  ■  C 


(21) 


This  equation  holds  for  the  forebody  and  afterbody. 

At  the  corner  of  the  body,  <5r  ■  0  but  6*  0,  which  means  the 
value  of  the  max-limm  thickness  Is  fixed,  but  Its  x  location  is  not 
fixed.  This  condition  gives 


A(F  -  r'Fr')  -  0 


(22) 


and  therefore 


<Co>  forebody  afterbody 


(23) 
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Ac  the  beginning  of  Che  forebody  (Che  subecrlpt  1  location  In 
Figure  2),  5x  -  0  end  6r  i  0.  Thla  glvee 


~Fr*  +  CPori  “  0 


(24) 


For  Che  forebody,  using  Modified  Newtonian  theory,  Equation  (2)  yields 


.*2 


S  "  CPo  1  +  r'2 


(25) 


and 


(Cpr')  -  C 


r,2(3  +  r,z) 
Po  (i  +  r»2)2 


(26) 


Now  use  equation  (21)  at  the  1  location  and  find* 


-r  (3  +  r  2) 

— - - - —  +  1 

(1  +  r^)2 


t 


CPo  /I  +  r'^z 


0 


(27) 


Using  a  Newton-Baphson  technique,  equation  (27)  can  be  solved  for  the 
optimum  initial  slope,  r^. 

Using  the  Modified  Newtonian  pressure  distribution,  it  follows 

that 


dCp 

dr' 


°Po 


(1  +  r'2)  2r’  -  r'2  2r’ 

(TTT^T2 - 


2CPo  r' 

"  (1  +  r'1^  ' 


(28) 


For  the  forebody,  equation  (21)  gives 


10 


t 

* 


I ,  SS  I'd 


r 


(29) 


2Cp  r,a  C£ 

(1  *  rU>'  ‘  C° 

or  rearranging 


t_  _  (1  +  r'2)2 _ 

Co  2C^r,  J  -  C£d  +  7rrPTZ  * 

Nov  on  the  forebody  section  x  can  be  related  co  r  by 


(30) 


(31) 

Applying  this  result  to  equation  (30) ,  It  follows  that 

1  dr  2Cpor'2  (1  +  r»2)  (r’2  -  3) 

Co  ***  [ICp^'3  -  Cf  (1  +  r’2)3^]  2 

Using  the  function  x/Cg  for  x  in  equation  (31)  it  follows  that 

2C  r’(l  +  r,2)(r’2  -  3)dr' 

,  x  po 

d  —  - - 

o  f 2C  r'  3  -  C.(l  +  r'2)  :,/2}2 

l  p0  £  J 

Equations  (30)  and  (33)  give  two  parametric  equations  to  determine 
the  optimum  forebody  shape. 

AFTERBODY 

The  afterbody  shape  optimization  starts  with  the  base  condition, 
fix  -  0,  but  fir  4  0  from  equations  (17)  and  (18) 


(32) 


(33) 


V 


3 

4  c  —  - 

2  PAB  rc 


0 


(34) 
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Differentiating  equation  (40)  by  r' 


dCa  PQ  d  |_P_ 

dr'  dr*  pQ 


pQ  d  f  p  1  dM2  d6 


q.  dM*  p  d6  drf 


Noting  that 


_ l 

cos^T 


\  +  i-=-±  M»]Y  ‘  1 

>  4 


v  JLm  2 

dC  P  Y  P 
p  o  _ _ _o _ 

de’  "  p-  ^FTT 


from  crvupsresiibie  aerodynamic  theory  subutition  into  equation  (4.1) 
yields 


dC 

P 

dC 

P 

d9 

V 

-£■  M2 

V  _  2Q 

dr' 

"  “die 

dr' 

"  "C 

CO  9  u  e 

rfT"=T 

Now  substituting  in  equation  (38)  the  optimal  equation  for  r  on  the 
afterbody  la  as  tallows 


P  Y  M2sin29 

r  _I_!o _ 

/M2  -  1 


- -  c 

r rnn  0 


Let  the  bracketed  quantity  in  equation  (43)  be  called  g(r')>  then 
C 

r  _ _ o 

r  g TFT 
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dr  o 


g'dr’ 


d  [JF-  - 

c„ 


Integrating  by  parts  from  the  corner  equation  (44)  becomes 


Chau 


dr,  .  m  -Ai1  -  1  d(M2) 

CO*“e  2M2(1  +  M2) 


and  cha  expression  for  x  on  the  afterbody  la 


x  -  x. 


C  C 
o  o 


r  t  .  rC|  _  I*1*  AP  -  \  d(M2) 

>  Co  C  t<2  Co  2sin20  K*(l  +  ^-r- 


o  +  *§*■  m2) 


The  Prendtl- Meyer  function  is  defined  a a 


8 


.  (47) 


■/ Y~J’  tan~ y'r^L-~p|M^  ~  ^  -  tan”1 45^  -  1  .  (48) 


Fur  comparlsion  with  Reference  5,  p  ■  should  be  true  at  0  -  0. 
The  Mach  number  when  0  «  0,  M^,  is  defined  by 


— .  (,  t-i 

P„  l1  +  2  **] 


or 


(49) 


*R 


2  . 


Y  -  1 


s 


-1 


with  po  being  the  stagnation  pressure  aft  of  the  normal  shock.  It 
should  be  noted  that  is  not  necessarily  the  same  as  Recalling 

equation  (38) 


8  -  -v  +  K 

o 

Ko  “  V(M-Ma) 


(38) 


v  ■  K  -  0 
o 


and 


for  6  lui  than  taro  on  tha  aftarbody. 

Tha  calculation  on  the  aftarbody  starts  at  the  corner  with  an 
assumed  8.  Equation  (45)  is  integrated  with  M  greater  than  until 

the  base  is  reached  where  equation  (37)  is  satisfied.  At  each  inte¬ 
gration  step  the  pressure  coefficient  is  defined  as 


C 

P 


p  p.. 

i 

o 


_____  _o 

p«  p- 


1 


(50) 


end  therefore 


C  * 


p  -  p 

ro 


1£  2 

q.  " 


or 


J_o 


C  +-2-1 


(51) 


Dividing  both  sides  of  equation  (51) 


(po/q-) 


-  I  + 


A  better  constant  can  be  defined  by  replacing  C  with 

o 


'i  c_ 


The  equations  for  the  forebody  using  Modified  Newtonian  pressure  theory 

are 
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On  the  afterbody  using  the  Prandtl-Meyer  expansion  for  tha  pressure 
coefficient,  the  equations  are 


c 


where 

9  -  K  -  v. 
o 

EULER IAN  RESULTS 

Equations  (52)  through  (55)  vara  digitized  in  a  marching  schema  to 
optimize  a  projectile  shape  given  Mach  number  and  initial  conditions. 

The  initial  slope  was  found  using  equation  (27)  which  started  the  marching 
schwa.  Transition  from  forebody  to  afterbody  was  made  at  the  maximum 
diameter  location.  The  end  location  (position  x ^  In  Figure  2)  was 
determined  by  minimization  of  base  drag  effects. 
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Casas  of  Mach  numbers  2  and  3  with  aea  lavel  conditions  wera  run 


and  compared  to  that  of  Reference  5.  Plots  of  thaaa  caaaa  ara  shown  in 
Figure  3.  Is  tha  Mach  3  case,  the  Eulerian  Optimization  gava  drag 
coaffidanta  which  wara  approximately  10X  less  than  tha  nethod  of 
Reference  5.  This  does  not  compare  wall  with  tha  experimental  data  of 
Reference  9  for  tha  Navy  25  mm  round.  Tha  wind  tunnel  data  in  Reference 
9  was  for  a  shape  optimized  on  a  L/3  ratio  of  5  with  maximum  diameter  at 
3.5  calibers  (x  /xe  -  .7)  at  Mach  3.  The  code  used  in  Reference  5 

C  £ 

predicted  a  drag  coefficient  22Z  less  than  experiment  while  the  Eulerian 
optimization  was  31Z  less.  The  Mach  2  case  is  much  worse  although  it 
should  be  noted  that  the  wind  tunnel  model  was  optimized  for  Mach  3.  A 
comparison  of  the  predicted  optimum  shape  for  L/D  «  3,  Mach  »  3  is  g±v«u 
is  Figure  4.  Both  shapes  are  very  similar  to  that  used  in  Reference  9. 

The  shapes  generated  by  Eulerian  optimization  were  similar  to  those 
generated  in  Reference  5  for  different  maximum  diameter  location.  The 
Eulerian  scheme  did  not  predict  the  drag  more  accurately  which  was  the 
goal.  Tals  led  one  to  question  the  accuracy  of  the  optimization  although 
the  shape  generated  was  essentially  the  same  as  previous  attempts.  The 
falling  of  the  drag  prediction  to  be  accurate  was  attributable  mostly  to 
the  wave  drag  prediction  of  the  modified  Newtonian  theory.  The  better 
accuracy  of  second-order  shock-expansion  was  then  given  consideration. 
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OPTIMIZATION  USING  SECOND-ORDER  SHOOT-EXPANSION  PREDICTION 
EULERIAN  TECHNIQUE 

The  second-order  shock-expansion  theory  was  used  by  Syvertson 
and  Dennis  In  Reference  10  to  predict  wave  drag  for  pointed  bodies  at 
angle-of-attack  equal  to  zero.  The  method  was  modified  by  Jackson  et. 
al.  In  Reference  11  to  account  for  blunted  bodies.  The  body  In  this 
method  is  replaced  by  a  tangent  body  which  is  a  series  of  conical 
frustrums  tangent  to  the  actual  body  at  various  body  locations.  An  attempt 
was  made  to  use  the  Eulerlan  optimization  scheme  with  the  second-order 
shock-expansion  method  as  the  wave  drag  component.  The  Newtonian  theory 
was  optimized  using  a  first-order  scheme  from  Reference  A.  A  second- 
order  scheme  also  from  Reference  4  was  initiated,  but  was  rejected 
due  to  complexity  of  the  terms  and  type  of  scheme  needed  for  numerical 
integration  of  the  shock-expansion  theory.  Of  major  concern  was  the 
large  size  of  required  matrix  operations.  A  less  complex  minimization 
scheme,  using  the  second-order  shock -expansion  method  for  the  surface 
pressures,  was  developed  and  is  described  below. 

SECOND-ORDER  SH0CR.-EXPANSI0N  PREDICTION  METHOD 

The  new  optimization  scheme  developed  here  is  essentially  a 
geometric  iteration  method  of  determining  an  optimum  chape.  The  accurate 
second-order  :>  Sack -expansion  technique  developed  In  Reference  17.  was 
chosen  to  replace  the  modified  Newtonian  plus  Prandtl-Meyer  expansion 
because  of  its  relatively  quick  computation  time  and  extensive  use  in 
body  alone  aerodynamics.  The  Mach  number  range  is  from  1.5  to  6.0  in 
this  method. 

The  original  second-order  shock-expansion  method  was  developed 
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for  pointed  noses  with  attached  shock  waves.  In  the  basic  method,  pressure 

on  the  Initial  cone  Is  obtained  from  a  cone  solution  and  Is  considered 

constant  on  the  cone.  The  pressure  drop  at  the  first  juncture  is  calculated 

from  standard  Prandtl-Meyer  expansion.  The  pressure  along  the  next 

frustrum  varies  exponentially  and  is  made  to  satisfy  three  boundary 

conditions.  The  first  boundary  condition  is  that  the  pressure  (p-)  just 

after  the  corner  of  the  initial  cone  and  first  conical  frustrum  Is  obtained 

from  Prandtl-Meyer  expansion.  The  second  boundary  condition  is  that 

the  pressure  gradient  (Sp/Ss^  at  this  position  (just  after  the  corner) 

is  obtained  from  an  approximate'  expression  developed  in  Reference  10. 

The  third  boundary  condition  is  defined  by  setting  the  pressure  at  infinity 

equal  to  the  cone  pressure  (p^)  that  would  exist  on  the  first  conical 

frustrum  if  it  were  infinitely  long.  The.  pressure  along  a  conical 

12 

frustrum  can  then  be  given  by 


P  -  pc  “  (PC  ‘  p2}  e 


where 


(56) 


[3p-]  (x  -  x.) 

—  - - — - — ■  •  — — -  (for  n  >  0)  (57) 

3*  J  2  <pc  "  P2J  COs52 

The  cone  angle  is  defined  aa  the  conical  frustrum  inclination.  T'-s 

Fre.ssure  gradient  just  downstream  of  the  comer  (position  2)  is 

10 

determined  from  the  approximate  expression 


where 


B 


ypM2 

2(M2-i) 


(59) 


X 


2  TP 
sin  2y 


n 


1 

M 


Jrfil 

_  2(y-l) 

1  +  M2 

“1ST — 

_  2  _ 


(60) 


(61) 


In  the  above  equations  (-3<$/3s)  is  the  curvature  of  the  surface 
which  is  zero  on  conical  frus trusts ,  ft  is  the  one  dimensional  area  ratio, 
and  the  subscript  1  refers  to  the  position  just  before  the  corner. 

Since  the  pressure  is  constant  on  the  initial  cone  (3p/3s)  equals  zero 
on  the  first  conical  frustrum  after  the  initial  cone.  For  all  sub¬ 
sequent  conical  frus tr urns  the  pressure  gradient  is  obtained  from  tne 
derivative  of  Equation  (56).  For  more  details  of  this  method,  see 
Reference  (10) . 

The  origiual second-order  shock-expansion  was  modified  by  Jackson, 

et.  al.**  for  blunt  bodies  by  using  the  modified  Newtonian  pressure 

distribution  up  to  a  "matching  point”.  Hie  matching  point  was  set  as  the 

maximum  angle  for  an  attached  shock  wave.  Beyond  the  matching  point , 

the  original  second-order  shock-expansion  is  used.  De  Janette  and 
12 

Jonea  made  two  modifications  to  that  of  Reference  11  that  increased 
accuracy.  A  computer  code  was  developed  using  these  modifications  along 
with  the  Van  Driest**  skin-friction  prediction  and  the  semi-empirical 
base  pressure  method  devised  by  Moore^. 


3S.  3!E^gSg35SHK3^JB!35SSSS^ 


Th*  modifications  Li&dts  in  Reference  12  consist  of  Introducing 

an  1  exact"  pressure  gradient  downstream  of  a  corner  and  determining 

a  new  matching  point  for  matching  second-order  shock-expansion  with 

modified  Newtonian  theory  on  blunt-nose  bodies.  The  "exact"  pressure 

gradient  is  baaed  on  the  method  of  characteristics.  The  following 

1 2 

equations  were  derived  on  the  surface  streamline 


Q  + 


U  sin  $ 


where 


Q 


~JiR  _ 

as 


(62) 


(63) 


Equation  (62)  is  integrated  numerically  around  the  corner  along  with 
the  Prandtl-Meyer  expansion  to  determine  Q.  Equation  C63)  is  solved 
for  the  pressure  gradient  (3p/3s).  A  new  matching  point  va9  found  to 
ba  the  position  on  the  nose  where  the  codified  Newtonian  pressure 
distribution  gives  a  local  Mach  number  cf  1.15. 


NEW  OPTIMIZATION  METHOD 

The  new  optimization  method  starts  with  a  semi-optimum  shape. 

An  iteration  method  is  then  used  to  determine  the  body  coordinates  which 
minimizes  the  total  drag  using  the  modified  second-order  shock -expansion 
method  to  calculate  surface  pressures  along  with  die  Van  Driest  skin- 
friction  and  empirical  base  drag  methods. 

The  selection  of  a  semi-optimum  body  began  with  a  review  of  the 


I 


4 


r*.  1 


optimum  studies  of  Reference  1  through  5  and  9.  The  general  conclusions 

of  the  first  five  references  Indicated  that  a  2/3  or  3/4  power  law  forebody 

gave  minimus  drag.  References  5  and  9  further  found  that  a  good  after  body  ..  . 

would  be  one  with  a  conical  boatail.  Further,  Reference  9  noted  that  for 

practical  applications,  a  blunted  nose  is  necessary.  A  review  of  experimental 

data  of  optimum  shapes  confirmed  the  theory  that  one  optimizes  for  a 

given  Mach  number.  An  arbitrary  selection  was  made  from  the  results 

of  the  above  study  that  the  semi “optimum  shape  would  be  chosen  for 

Mach  3.  The  initial  bluntness  was  made  L/D  dependent  from  the  results 

of  the  Euleriaa  Optimization.  The  forebody  was  set  as  a  3/4  power  law 

body  allowing  for  the  selection  of  different  maximum  diameter  positions. 

From  the  maximum  diameter  location  aft  a  6°  conic  was  chosen.  The 

resulting  semi-optimum  body  differed  from  Reference  5  and  the  Eulerlan 

Optimization  in  che  forebody  shape  and  the  boatail  cutoff  location. 

A  total  of  20  coordinates  were  selected  as  an  adequate  description  of  the 

body  with  14  on  the  forebody  and  the  remaining  6  on  the  boatail. 

20 

A  set  of  independent  coordinates  {  x.  }  were  selected  with  the 

1  i-1 

i»14  point  taken  as  the  point  of  maximum  diameter  (note  i  in  this  section 

represents  a  coordinate) .  The  corresponding  set  of  dependent  variables 
20 

{  r.  }  were  initially  determined  from  the  semi-optimum  body.  For 

1  i-1 

prescribed  values  of  {x^},  r^»  and  Reynolds  number,  the  drag 
coefficient  can  be  represented  by 

CD  -  (^(r^  i  -  1,2,..., 20  (1*14) 

It  is  desired  to  determine  the  values  of  r^  which  makes  a  relative 
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The  derivatives  in  Equation  (67)  are  formed  by  the  following  central 
difference  quotients 


(68) 


(69) 


where 


S*  '  Vrj,o’  ri.o  +  Sri> 

\  m  Vrj.° •  ri.o  '  Srl> 


1,2 . 20;  j  i  i 


and  ''r  is  two  percent  of  r^. 

The  iteration  process  involves  calculating  CQ  ,  C+  and  C~ 

o  i  i 

1,..,,20).  Then  using  equations  (68)  and  (69),  equation  (67)  can 
b*  sed  to  calculate  Ar^.  Then  new  values  of  r^  are  calculated  by 
adding  Ar^  to  the  old  values.  The  iteration  process  is  continued 
iint.il  '.onvergence  which  was  assumed  to  occur  when  Ar^  changed  less  than 
one  ccent.  In  the  iteration  process,  if  |Ar^|  >  5r^  then  the  magni¬ 
tude  of  Ar^  was  taken  to  be  Ar^.  Again,  note  that  r^  for  i«14  was  not 
changed  since  it  is  the  maximum  diameter  point. 

Convergence  did  not  occur  in  cases  where  the  maxiimim  diameter 
location  was  less  than  25Z  of  the  total  length.  This  is  probably  due 
to  the  negligence  of  the  cross  product  terms. 
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RESULTS 


The  optimization  scheme  described  above  was  digitized  in  an 

efficient  maimer  to  minimize  computation  time.  Since  the  maximum 

diameter  location  is  an  input  parameter,  cases  vere  run  varying  this 

location  (x  lx,).  Also  varied  were  Mach  number  and  total  length  to 
c  t 

diameter  ratios.  Results  of  typical  runs  for  sea  level  conditions 
are  given  in  Figures  5,  6,  and  7.  As  a  comparison,  the  shape  generated 
by  the  new  optimization  technique  is  drawn  with  that  of  Reference  5  in 
Figure  8.  This  shape  also  is  very  close  to  the  shape  generated  by  the 
Eulerlan  technique  (Figure  4).  The  predicted  drag,  however,  is 
different.  The  Mach  3,  L/D  equal  5  case  is  found  in  Figure  9.  The 
Mach  3  Experiment  point  is  that  found  in  the  wind  tunnel  tejt  of 
Reference  9.  The  25  mm  shape  tested  is  quite  similar  to  those  in 
Figure  8  with  the  exception  of  grooves  placed  on  the  boatail.  These 
grooves  are  used  for  rotating  band  attachment  and  could  be  responsible 
for  some  of  the  9.6Z  difference  in  drag  coefficient.  The  shapes  for 
other  cases  using  the  new  scheme  compared  similarly  for  other  Mach 
numbers,  that  is  good  shape  agreement,  but  different  drag  coefficients. 
An  interesting  development  in  this  method  is  that  the  design  curves 
produced  are  flatter  in  the  optimum  dreg  area  than  those  of  Reference  5 
This  would  tend  to  give  projectile  designers  more  freedom  in  actual 
shape  variation  and  still  produce  low  drag  results.  A  summary  of  the 
Mach  3  cases  are  shown  in  Figure  10.  These  curves  Indicate  the  trend 
of  Increased  x  lx,  with  decrease  of  L/D  for  optimum  drag.  Figures  5, 

C  £ 

6,  and  7  illustrate  the  trend  of  increase  in  x  lx.  with  Increase  in 
design  Mach  number. 

The  new  optimization  Iteration  code  is  simple  to  operate  and  gives 
the  user  ease  in  running  multiple  cases.  The  number  of  iterations  to 
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convergence  ranged  from  A  to  17 .  An  average  case  (1  Mach  number  and  L/D) 
coat  approximately  $10  on  both  the  IBM  370  and  CDC  6700.  Core  require¬ 
ments  are  minimal  and  the  code  could  be  put  on  larger  mini-computers 
(64K  bytes).  Output  includes  the  number  of  iterations  to  convergence , 
components  of  drag,  total  drag  coefficient,  and  the  minimum  drag  shape 


coordinates. 


CONCLUSIONS  AND  RECOMMENDATIONS 


1.  Two  numerical  methods  vara  davalopad  for  calculating  optimal 
projectile  shapa  for  gtin-tmtn*  total  drag. 

2.  Tha  Eulerl&n  optimization  technique  calculates  similar  shapes 
for  minimum  drag,  but  is  inaccurate  in  its  prediction  of  total  drag. 

3.  Tha  new  optimization  technique  gives  both  an  optimum  shapa 
and  a  more  accurate  drag  prediction  whan  compared  to  experiment. 

4.  A  limitation  of  the  new  optimization  code  la  that  the  maximum 
diameter  location  must  be  greater  than  25X  of  the  total  length. 

5.  The  ratio  of  maximum  diameter  location  to  total  length  tends 
to  Increase  for  decreasing  L/D  ratio  and  Increase  with  Increasing  Mach 
number  for  optimum  shapes. 

6.  A  good  agreement  between  three  different  predictive  techniques 
lends  credibility  to  the  actual  shape  of  minimum  drag  rounds. 

7.  This  technology  should  be  proved  experimentally  in  both  large 
caliber  rounds  (such  as  76  mn)  as  well  as  small  caliber  rounds. 
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Figure 


Comparison  of  Drag  Coefficient  va  L/D  for 


EULERIAN  OPTIMZATION 


FIGURE 4.  OPTIMUM  SHAPE  COMPARISON, 

EULERIAN  OPTIMIZATION  AND  REFERENCE  5. 
MACH  NO.  =  3  ,  ^  =  5 


Coefficient  vs  x  /x.  for  L/D-6,  New  Optimization 


New  Optimization 


A  L/D-5 


Figure  10.  Drag  Coefficient  vs  x  /x  ,  Mach  No. *3,  IJew  Optimization 
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drag  coefficient; 

base  drag  coefficient 

akin  friction  coefficient 

mean  akin  friction  coefficient 

pressure  coefficient 

base  pressure  coefficient 

diameter 

reference  diameter 

length  of  configuration 

Mach  number 

pressure 

Frandtl  number 

dynamic  pressure 

radius  of  body  (Figure  2) 

body  slope,  dr/dx 

Reynolds  number 

turbulent  recovery  factor 

reference  area 

vetted  surface  area 

temperature 

wall  temperature 

velocity 

length  coordinate  (Figure  2) 

ratio  of  specific  heats 

angle  along  body  surface  (tan_1(dr/dx)) 


GLOSSARY  (Cont'd) 


y  coefficient  of  absoluce  viscosity 

v  Prandtl “Meyer  function 

0  density  of  air 

n  ratio  of  cross-sectional  area  of  streamtube  to  that  at  M*1 

Subscripts 

0  stagnation  condition 

1  condition  immediately  before  a  corner 

2  condition  lnsnediately  after  a  corner 

B  base  conditions 

c  position  of  maximum  c  '.ameter  (Figure  2) 

f  position  at  end  of  body  (Figure  2) 

1  position  at  front  of  body  (Figure  2) 

«  freestream  conditions 
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